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Abstract 

We give a theoretical treatment of magnetic dipole absorption of electromagnetic 
radiation in small conducting particles, at photon energies which are large compared to 
the single particle level spacing, and small compared to the plasma frequency. We discuss 
both diffusive and ballistic electron dynamics for particles of arbitrary shape. 

The conductivity becomes non-local when the frequency is smaller than the frequency 
uj c characterising the transit of electrons from one side of the particle to the other, but in 
the diffusive case u> c plays no role in determining the absorption coefficient. In the ballistic 
case, the absorption coefficient is proportional to u 2 for u> <C u> c , but is a decreasing 
function of u for u ^> u c . 



1 Permanent address: Department of Physics and Applied Physics, John Anderson 
Building, University of Strathclyde, Glasgow, G4 ONG, Scotland, U.K. 

2 Permanent address: Max-Planck-Institut fur Physik komplexer Systeme, Nothnitzer 
Str. 38, 01187 Dresden, Germany 



1. Introduction 



For sufficiently low frequencies, the interaction of small conducting particles with 
electromagnetic radiation is dominated by absorption rather than scattering. The classical 
theory of the interaction of electromagnetic radiation with spherical particles of uniform 
composition was considered by Mie [1]; the theory encompasses conducting particles with 
a complex dielectric constant, which is often modelled by the Drude theory [2]. If the 
particles are small compared to both the wavelength and the electromagnetic skin depth 
of the radiation, the dominant contributions to the absorption are called the electric and 
magnetic dipole terms [3]. In both cases the absorption is due to Joule or Ohmic heating 
caused by electrical currents flowing through the particle: the electric dipole term is due 
to currents which establish electrical polarisation of the particle, and the magnetic dipole 
term is due to eddy currents induced by variation of the magnetic field. 

At frequencies below the plasma frequency, the electric field is screened from the 
interior of the particle, but the magnetic field can penetrate the whole of the particle. 
Although magnetic effects are negligible in atomic absorption processes, they could become 
significant when the number of atoms in the particle sufficiently large that most of the 
atoms are screened from the electric field. In fact, magnetic dipole absorption is often 
the dominant absorption process in suspensions of small metal particles [4]. Very few of 
the many theoretical papers on absorption of radiation by small particles have considered 
the magnetic dipole contribution, some exceptions are [5,6] which consider magnetic dipole 
absorption in the context of effective medium theories. Because it is typically the dominant 
contribution, it is appropriate to consider the problem of magnetic dipole absorption in 
some detail. 

The Mie theory is restricted to spherical particles in which the dynamics of the charge 
carriers is diffusive, and quantum mechanical effects are ignored. In this paper we will 
give the first theoretical treatment of magnetic dipole absorption going beyond the Mie 
theory: we will consider both diffusive electron dynamics (for which the particles are char- 
acterised by their bulk conductivity), and ballistic electron motion (in which the particle 
is smaller than the bulk mean free path of the conduction electrons) . Our approach allows 
for arbitrary particle geometries, and we give careful consideration to the fact that the 
conductivity is non-local when the particles are very small: quantum mechanical effects 
are included using a semiclassical approach. The paper complements [7-9], which gave a 
comparably comprehensive treatment of electric dipole absorption. 

We do not explicitly consider the structures in the absorption close to the single 
particle level spacing which were originally considered by Gorkov and Eliashberg [10]. 
These structures are determined by repulsion between energy levels and the appropriate 
tool to analyse them is random matrix theory. Their full characterisation requires an 
estimate of a mean-square matrix element, which was not given correctly in [10]. Our 
results for the low-frequency limit provide the correct estimate of this quantity for magnetic 
dipole absorption. 

Sections 2 to 4 will be concerned with various aspects of the formulation of the prob- 
lem, discussing respectively the relation between the absorption of radiation and correlation 
functions of the electron motion, the criteria for a self-consistent solution of the equations 
determining the electric field driving the eddy currents, and the definition and semiclassical 
estimation of the non-local conductivity which is required to determine the self-consistent 
field. Our semiclassical estimates for the non-local conductivity are closely related to ex- 
pressions given by Argaman [11]. In section 5 we develop the theory for magnetic dipole 
absorption in particles with diffusive electron motion, assuming that the electric field is 
known. This calculation also yields the form of the non-local conductivity applicable to the 



diffusive case: in section 6 we consider the solution for the self-consistent field, and discuss 
results for some specific geometries. Our formula for the non-local conductivity is identical 
to one given by Serota and co-workers [12,13], who used diagrammatic techniques. Our 
derivation is more direct and requires fewer assumptions: we discuss this point further in 
section 6. 

Our approach uses a semiclassical estimate described in [14], which relates mean 
squared matrix elements to classical correlation functions. We might expect that there 
should be features in the absorption spectrum which are related to the characteristic 
timescale for decay of classical correlations, in this case the typical time for a particle 
to cross the particle: the importance of this timescale was emphasised by Thouless [15], 
and in the case of diffusive electron motion, we will refer to the characteristic frequency 
scale uj c = D/a 2 as the Thouless frequency (D is the diffusion constant and a is the charac- 
teristic size of the particle) . Another reason for expecting uj c to play a role in determining 
the absorption coefficient is that the conductivity is non-local when ui is not large compared 
to uj c . We find however that u c plays no role in the final expression for the absorption 
coefficient. 

We consider the case of ballistic electron motion in section 7. We are only able to gain 
limited information about this case: we find that the absorption coefficient is proportional 
to uj 2 at frequencies small compared to u c = vf/ci, and that it is a decreasing function of 
frequency for u> ^> u c . 

Our conclusions are in many ways parallel to those for electric dipole absorption. In 
ballistic systems, it was found [7,8] that the electric dipole absorption has resonances in 
the absorption coefficient with a frequency scale u c = a/vp, where v-p is the velocity at the 
Fermi energy. By contrast, in the case of diffusive electron motion, it was found [9] that 
there is no structure in the absorption coefficient at the frequency scale uj c = D/a 2 . 

Finally, we remark that there is a large literature concerned with the effects of time 
dependent magnetic fluxes on metallic loops: when the magnetic flux varies sinusoidally, 
the absorption of energy by the loop is a special case of the magnetic dipole absorption 
which we consider here. Most of the papers on this topic are concerned with quantum size 
effects analogous to those considered by Gorkov and Eliashberg [10]: two recent works in 
this area are [16] and [17]. 



2. Formulation of the problem 



The absorption of radiation is usually described by an extinction coefficient 7(0;), 
which is defined as the fractional loss of intensity per unit length of sample, divided by 
the volume fraction F occupied by the particles. We will express our results in terms of 
the rate of absorption of energy (dE/dt) within a single particle. If the amplitude of the 
electric and magnetic fields are Eq are Bq respectively, the intensity of the radiation is 
I = ^qEq = and the relationship between 7 and (dE/dt) is therefore 

2 /dE\ . , 

where V is the volume of a single particle. In this paper we will define the absorption 
coefficient a(ui) as the rate of absorption of energy for a single particle, divided by the 
electric field intensity: 

, N 1 /dE\ 1 /dE\ , nn . 

The normalisation with respect to electric (rather than magnetic) field intensity is used to 
facilitate comparison with the results in [7-9]. 

The particle will be considered to consist of a static potential well which traps a gas 
of non-interacting fermions (electrons), initially with occupation probability f(E) (which 
would be identified with the Fermi-Dirac distribution). In a quantum mechanical calcula- 
tion the rate of absorption of energy is determined by the Fermi golden rule. This states 
that the rate of transition under the action of a periodic perturbation with frequency uj 
and matrix elements AH nm , from a state with energy E n , to a quasi-continuum of final 
states with energies close to E m = E n + Tiu> is 

R=^g(E m )(\AH nm \ 2 ) u , (2.3) 

where g(E m ) is the density of final states and ((Ai^™! 2 )^ is the mean-square matrix 
element for transitions from E n to states close to E m . The energy absorbed by an electron 
making an upward transition is Tiu). The rate of absorption of energy is therefore 

(^j = huj J dEg(E) R(E) [f(E) - f(E - hu)] ~ RghW (2.4) 

where the approximate equality is applicable in the limit where fuo and kT are both large 
compared to the mean level spacing, but small compared to other energy scales: in the 
right hand expression both g and R are evaluated at the Fermi energy. 

The mean-square matrix element can be estimated semiclassically [14] : 

1 r°° 

(|Atf nm | 2 > w = — — / dt exp(iut)(AH(t)AH(0)) E (2.5) 

where (AH(t)AH(0)) e is the micro canonical autocorrelation function of the classical ob- 
servable corresponding to AH, evaluated at energy E. Combining this result with (2.4), 
the absorption coefficient can be expressed in terms of the classical autocorrelation func- 
tion of the perturbation Aff (r, p). The resulting expression can also be written in terms 



of the classical change in energy of the individual electrons due to the perturbation: the 
change in energy of an electron following a trajectory r(£), p(t) is 



AE(t) = J dt' — (v(t'),p(t')) . (2.6) 

Combining this result with (2.5), the rate of absorption of energy by the electron gas is 

dE\ g d 



where (AE(t) 2 ) is the variance of the change in the single-electron energies. 

In our problem the perturbation is a sinusoidally varying electromagnetic field, spec- 
ified by a vector potential A(r) exp(iut), and a scalar potential <fr(r) exp(iut). The com- 
ponent of the perturbation of the Hamiltonian which is quadratic in A can be neglected 
when calculating the leading order absorption coefficient; the remaining terms are 

AH= -iL(p.A + A.p) +e$ . (2.8) 

In discussing the magnetic dipole absorption, it is given that the magnetic field is 

B(t) = V A A = B exp[iut] . (2.9) 

The fluctuating magnetic field induces an electric field E, which is given by 

dA 

E=-— + V$. (2.10) 
ot 

The Hamiltonian admits a set of gauge transformations (A, $) — > (A', $') = (A, $) + 
( V/U, dtn) which leave the electric and magnetic fields unchanged. We will assume that the 
gauge has been chosen so that $ = 0. Physically, the electric field is not uniquely defined 
by the magnetic field, and must be determined by a self-consistent condition, which we 
discuss in section 3. 



3. Self-consistent choice of the vector potential 



3. 1 Self-consistent electric field 

The eddy currents induced by the fluctuating magnetic will themselves generate a 
magnetic field. We will consider only the case of extremely small particles, for which this 
additional magnetic field is negligible. We therefore assume that the magnetic field is 
simply the externally applied field, B(£) = Bo exp[iu;t]e3. The electric field is required to 
uniquely determine the perturbation of the Hamiltonian. In this section we will discuss 
the self-consistent calculation of the electric field. 

In discussions of the Zeeman effect in atomic physics, the perturbation of the Hamil- 
tonian representing the magnetic field is conventionally taken to be proportional to the 
component of the angular momentum operator along the direction of the field. It is natu- 
ral to ask why a more involved procedure is used here, but is unnecessary for the Zeeman 
problem or for calculation of static magnetic susceptibility. At the end of this section we 
show that the angular momentum operator gives the correct answer in the limit where the 
frequency approaches zero, but not in general. 

The electric field satisfies the Maxwell equations 

VAE=^ V.E=^. (3.1) 
at e 

The electric field causes a current density j to flow within the particle. We will assume that 
a linear response theory is valid, but in general the current may be a non-local function of 
the electric field: we will write 

j(r,t)=fdr'f dt'a(r,r';t-t')B(r',t') (3.2) 

J J — oo 

where a is the non-local conductance tensor. We will also write this relation in the form 

j(r,u;) = J dr' a(r,r';u)E(r',u>) = aE(r,u) (3.3) 

where the second equality defines an operator a which maps the electric field E(r, ui) non- 
locally into the current field j(r, uS). For a monochromatic perturbation, the charge density 
is 

P = --V.j. (3.4) 

Combining these results, we find the following equation for the electric field 

V.j' = 0, j' = (a - iu;eo)E . (3.5) 

This equation must be supplemented by a boundary condition in order to uniquely deter- 
mine the electric field. This is 

n.j = (3.6) 

where h is a unit vector normal to the boundary: this condition represents the fact that 
charges cannot enter or leave the sample. 

3.2 Representation of the field in terms of potentials 

It will be convenient to write the electric field in terms of a vector and a scalar 
potential: 

E = VAV + V(^>. (3.7) 



We will only consider in detail cases where the field ip(r) is of the form 

V = V(^2/)e 3 - (3.8) 

This form is appropriate when the conducting particle is two dimensional, lying in the 
plane z = 0, for three dimensional particles in the form of general cylinders aligned with 
the z axis, and can be extended to spheres and some other non-cylindrical geometries. 
Substituting (3.7), (3.8) into the Maxwell equations, we find that ifi satisfies Poissons's 
equation in the form 

V 2 V> = iujB . (3.9) 

We will always choose ip(x, y) to satisfy the condition that ^ = 0on the boundary. Hav- 
ing uniquely specified ^(x,y), equations (3.5) and (3.6) are transformed into a equations 
determining the scalar potential (p. 

In some cases a local, isotropic conductivity S(u;) will provide an adequate description. 
In this case, the condition (3.5) reduces to the requirement that V.E = 0, in which case 
the electric field can be written in the form 

E = VA<0 + V(^>, V 2 (/> = 0. (3.10) 

The boundary condition corresponding to (3.6) is then satisfied by taking a solution for 
which ip = and <j> = on the boundary: the latter condition implies that 0=0 every- 
where. 

3.3 A remark on the low-frequency limit 

After having described the approach used to define the correct perturbation, we will 
now show that any form of the electric field which has a uniform value of V A E gives 
the correct value of the absorption coefficient in the limit uj — > 0. According to (2.7), the 
absorption coefficient is proportional to the variance of the change in the single particle 
energy. The change of the single-particle energy can be written 

AE(t) = [dr.E= fdt' ^j. E(r(£')) = fdt' ^. A(r(t')) • (3.11) 
J Jt J to dt 

Consider the effect of making a transformation A — > A' = A + Wip on the absorption 
coefficient. The change in the single particle energy is transformed to AE': 

ft dr 
AE' = AE + iuj / dt' exp(iujt') — .V(p 

= AE + iu [ dip[v(t')] exp(iwt') = AE + iuX(t, u) . (3.12) 

Jto 

If ui = 0, the correction X(t) introduced by the transformation is simply 

X(t, 0) = jf dt' —.V<p = <p[r(t)] - ^r(to)] (3.13) 

which remains bounded as t — > oo. For finite u the correction X(t, uj) satisfies 

(X\t^)) = t J°° dr expM^(r)^(0)^+O(l)=rHi + O(l) (3.14) 

provided the correlation function of d<p/dt decays faster than r _1 . Comparison with (3.13) 
shows that the coefficient r(o;) approaches zero as uj — > 0, implying that the gauge depen- 
dent contribution to the absorption coefficient vanishes in the limit ui — > 0. 



4. Semiclassical theory for non-local conductivity 



4.1 General formula 

We will use a semiclassical analysis for the non-local conductivity <r(r, r',u>). We 
will first consider the problem in rather abstract terms: we will discuss a Hamiltonian 
H(r, p,X), where X is a time-dependent parameter. The Hamiltonian determines the 
motion of particles in a gas with phase space density p(r, p, t). The phase space density 
satisfies the Liouville equation dtp = {p, H}. A solution can be written in the form 

p(r,p;t) = f(H(r,p;X)-E F ) 

-X fjt' ff(r(0,p(0;*)f^(#(r,P;* - E F ) +0(X 2 ) . (4.1) 

Formally, this is an expansion in the velocity of the perturbation, X: the results will be 
valid for all frequncies, because the amplitude of the perturbation is infinitesimal. 

The leading order Weyl or Thomas-Fermi estimate of the density of states of a quan- 
tum system states that the density of quantum states is (27rh)~ d in classically accessible 
regions of phase space. We will therefore multiply the above solution by this factor, and 
take f(E) to be the Fermi-Dirac function. 

Our Hamiltonian, H = (p — eA) 2 /2m e + V has a time dependent vector potential, so 
that we can write 

• dH dH „ e „ 
X— = —.B= — p.E. (4.2) 
OX OA m e V ; 

The resulting current is 

j(r,£) = — f dppp(r,p,*). (4.3) 
The current flowing in response to the electric field E(r, t) is therefore 



e 2 



{2ith) d m, 



fdt'J dp J dr / ^p i P i (r,p;t , -t)*[r , -R(r,p;t-0]S i (r / ,0 (4.4) 



where Pj(r, p; r), i?i(r, p, r) are the i th components of the momentum and position at time 
r for a trajectory which starts at (r, p) at time t = 0. The components of the non-local 
conductivity tensor are therefore 

Vij(r,r';t) = (2 J^ m 2 / ^P / dr'^p 4 P i (r,p;0*[r'-R(r,p;t-0] • (4-5) 
We can write this result in a simpler form 



where the angle brackets denote an average over the initial momenta, defined by (4.5). 
We will consider the evaluation of this quantity for diffusive motion in section 6; next we 
consider the case of ballistic motion. 



4.2 Results specific to ballistic systems 



Equation (4.5) can also be expressed as a sum over classical trajectories which travel 
between r and r': 



paths 



det 



dR k 
dpi 



(Pinit)j(Pfin)i • 



(4.7) 



In the low temperature limit the term df/dE reduces to a delta function, and this expres- 
sion becomes a sum of delta functions 5(t — Tj), where the Tj are the times of trajectories 
from r to r' at the Fermi energy Ep. 

It is more convenient to consider the frequency dependent non-local conductivity: if 
the electric field is E(r) exp(io;t), then the current can be written in the form 



ji(r,t) = exp(iut)^2 / dr'aij(r,r';u) E 3 (v') 

i 

Comparing with (4.4) and (4.5), we find 

^(r, r » = /drexp M o,,,,- , 



(4.8) 



j dr exp(i^r) J d6 Pi (0) p'^9) 8(r' - R(r, p; r)) (4.9) 



(2nK) 2 m 2 e 

x ' e paths 

where in the second line we have specialised to the case of two dimensions, and 6 is the 
initial angle of the trajectory. Performing the integrations, we find 



0-y(r,r';u;) = 



E 



Pi(Q)p'j(Q) exp(iwr fc ) 



(4.10) 



where the sum runs over all trajectories which travel from r to r' in time tu at the Fermi 
energy. This general expression can be specialised in a variety of ways. We remark that it 
has a rather simple form for billiards with boundaries consisting of only straight edges. In 
this case the times are proportional to the lengths of the trajectories, and because 
there is no focusing or de-focusing of bundles of trajectories when they bounce off the 
boundary, the form of the determinant is very simple: we find 



<7ij(r, r';u) 



Pf^L^, 1 nf'ny exp(iu)L k m e /p F ) 



(2nTi) 2 m 



(4.11) 



where n^, n'j are the components of a unit vector in respectively the initial and final 
directions of the trajectory. As an example, the case of a square billiard is illustrated in 
figure 1. The formula also gives the non-local conductivity in free space, with only the 
direct trajectory included. 

We close this section by remarking that in the limit u> — > oo the non-local conductivity 
becomes a very rapidly varying function of r', except for when the path length of the 
trajectory is very short. Unless the electric field is a rapidly varying function of position, 
the dominant contribution to (4.4) comes from the region where r' is close to r, implying 



that a local conductivity £y (r, u>) will give an adequate description. When r' is close to r, 
the non-local conductivity can be approximated by (4.11), with only the direct trajectory 
included. The local conductivity is then obtained as follows: 



(? f 1 



ie 2 Pp x iNe2 x (a -\o\ 

— = ~ ( 4 - 12 ) 



2nh m uj m e u 



where N is the electron density per unit volume. This result is precisely the same as the 
high frequency limit of the Drude formula for the conductivity. 



5. Absorption coefficient for diffusive electrons 



5.1 Preliminary comments 

In the present section we calculate the absorption coefficient ot(u>) for systems with 
diffusive electron motion, using (2.6) and (2.7), assuming that the self-consistent electric 
field is known. We show that the absorption coefficient can be written as a sum of two 
terms. The first term describes a classical bulk contribution. The second term introduces 
boundary contributions which could modify the absorption coefficient at frequencies below 
uj c = D/a 2 . 

We begin by briefly discussing the classical expression for the absorption coefficient: 
it is natural to compare the final answer with this result. The rate of absorption of energy 
is given by integrating the rate of Joule heating j.E over the volume of the particle. The 
current density j is proportional to the local electric field: j = EoE = kuEoA where Eo is 
the bulk conductivity of the metal. We therefore have 

v h |j|2 ^ h A2 ^ w / dt A2 (8 - 1) 

where n is the density of states per unit volume, n = dN/dE = g/V. The boundary 
condition for the electric field is determined by the fact that the current must be tangential 
to the boundary. Unless there is a constant biasing magnetic field present, E is aligned 
with j. This implies that E.n = at all points on the boundary (where n is a unit vector 
normal to the boundary). 

Equation (2.5) suggests that the absorption coefficient might exhibit deviations from 
classical behaviour at frequencies small compared to the Thouless frequency uj c = a 2 /D : 
which is the inverse of the time taken for an electron to diffuse across the sample. In 
the following, we give a semiclassical treatment of the absorption coefficient with diffusive 
electron motion. There are two distinct but related issues which must be addressed here. 
Firstly, for a given electric field E(r), does the absorption coefficient exhibit any structures 
at the Thouless energy? Secondly, is the self consistent solution for the electric field 
different above and below the Thouless frequency? In this section we consider the first 
of these issues. In section 6 we will show how one of the results below can be used to 
determine the non-local conductance, and consider the determination of the self consistent 
field in greater detail. 

5.2 Calculation of the energy absorbed 

Our calculation is based upon (2.6) and (2.7). Because the instantaneous velocity is 
not well defined for a diffusive trajectory, we will divide the trajectory of the electron into 
finite segments, in which the electron travels from r n to r n+ i with a uniform velocity, in 
a fixed time increment dt. The r n are chosen from an ensemble of random walks confined 
within the boundary of the particle. The change in the single-electron energy is 



AE(t) = Re 



Re 



leu 



[ l dr' 
/ dt' ex.p(iut)A(r).— 
Jo dt 



A/--1 



(5.2) 



ieu ^2 exp(icut n )A(r n ).5r n 

where t = Afdt, t n = (n + ^)5t, 5r n = r n+ i — r n , and the quantity A is defined by the 
requirement that each term in the sum equals the contribution to the integral from the 



corresponding link in the random walk. For diffusive motion with fixed diffusion constant, 
we must take (Sr 2 ) ~ St, so that the error in each step must be 0(Sr 3 ) to achieve a 
convergent estimate of the integral. The sum (5.2) can be approximated as follows: 



AE(t) = Re 



ieuj ^2 exp(iwt n )A[i(r n + r n+ i)].(r n+ i - r n ) 



n=0 



+ 0(AfSr 3 ) (5.3) 



and the error term vanishes in the limit St — > 0. We remark that equation (5.2) is a 
stochastic integral, and evaluation of the integrand of the midpoint of the step is equivalent 
to using the Stratonovich definition of the integral [18]. The variance of (5.3) is 

Af-l 

(AE(t) 2 ) = eV J2 exp[ia;(t i -t J -)]<A[i(r i + r i+ i)].(5r i A[i(r J - + r j+1 )].S rj ) + 0(St) 

i,j=0 

oo 

~ e 2 u 2 M ^2 ex.p(iut n )C n (5.4) 

n= — oo 

where 

C n = (A[±(r + ri)].5r A[|(r n + r n+1 )].Sr n ) (5.5) 

The absorption coefficient is therefore proportional to the Fourier transform of a correlation 
function C(t): 

a ^ = 25tB^ exp(iwt n )C n = -^r dtexp(iut)C(t) (5.6) 

^= — 00 J — oo 

where C(nSt) = C n /(St) 2 . The correlation function can be expressed in terms of the 
propagator P(r, r'; t) which gives the probability density for reaching r' from initial position 
r at time t: for times large compared to 8t this satisfies the diffusion equation [dt— DW 2 ]V = 
0. Assuming summation over repeated indices, we find 

C{t) = ^ / dr J dr ' J d6r J d8r'Ai(r+ ±5r) Aj(r' + \8v') 6n Sr'j 

xV(r + St, r'; t - St) V{r, r + Sr; St) V{r', r' + Sr'; St) . (5.7) 

We discuss the cases t = and t ^ seperately. At t = 0, we have Sr = Sr' , and the 
correlation function is easily evaluated, giving a result which is of order 0(St): 

C = ((A.Sr) 2 ) = ^ J drA(r) 2 . (5.8) 

Here we have used ((A.Sr) 2 ) ~ ^(A 2 )(Sr 2 ) = DSt(A 2 ). This result implies that there is a 
contribution to C(t) of the form (D/V)S(t). 

The case t ^ requires a more delicate treatment. We can expand (5.7) in both the 
short time interval St and in the short step Sr. It will turn out that the leading order 
contribution is of the order of 0(St 2 ), as opposed to the case t = 0. Because the motion is 
diffusive, we have Sr 2 ~ St, so that terms up to quartic in Sr must be retained 

C(t) ~ ^ J dr J dr' J dSr J dSr' [A t + \d rk A % Sr k + ^ ri A % Sr k S n ] 



x [A'j + \d r , k A] 8r' k + \d\ r[ A^ 8r' k 8r[] 8 n 8r' 3 

x [P(r, r'; t) - d t V(r, r'; t) 8t + d rk V(r, r'; t) 8r k + |<9 r 2 fe7 ,P(r, r'; t)8r k 8 n ] 

xP(r, r + 8r; St) V(r, r' + 8r; 8t) . (5.9) 

The terms containing d t V(r, r'; t) can be dropped when there are more than two factors 
of Sr. The integrals over products of the 5r can now be separated out to give 

C{t) = J^JdrJ dr'^A i A j [P-dtVSt}{8r i ){8r' j ) 

+[\A' j d rk A i V+\A' j d rk {VA i ) - lA'^AiV^SnSrk)^) 
+[±A i d rk A>V](8r' j 8r' k )(8r i ) 
+ [±Ad 2 rkr A' 3 V + \A % A'^VKSriSrkSn)^) 

+[M 1 ^ r A;?](fr;Kfr;)(fr ! ) 

+[|a n ^a rfc ^P+ lAid^AjdrrfiSnSrkXSrlSr^ (5.10) 

where (8r) = J d5r 5r P(r,r + 5r; 5t). Now consider the form of these integrals when 8t is 
sufficiently small. The propagator "P(r, r + 8r; 8t) is small unless 8r is small. When r is not 
close to the boundary, this propagator can be approximated by a function of the distance 
travelled, Po(|5r|, St): because the steps are assumed to be independent, the variance (8r 2 ) 
averaged over this distribution can be identified with 2dD8t (where d is the dimensionality 
of space). When r is close to the boundary, a solution satisfying the boundary condition 
n.VP = is constructed by the method of images. We denote the image of the source 
point r by r* . The diffusion propagator is then 

V(r,r + 6r;6t)~Vo(\6r\,6t) + V (\6r + T-r*\,6t) . (5.11) 

Since 8t is small compared to the time to traverse the particle, the second term only 
contributes for points r close to the surface, which can thus be considered locally flat. 
We then introduce a local coordinate system arranged so that the nearest boundary point 
defines the origin. In two dimensions, the surface tangent is given by the line x = 0, and 
the normal by y = 0. the point r lies at (x, 0), and r* = (— x, 0) (figure 2). 

The average (Sri) vanishes unless r is close to the boundary, in which case the mean 
displacement is inwards, and its projection in the direction perpendicular to the surface is 

P oo 

(8x) x = dx'(x' -x)[f(x' -x) + f(x' + x)} (5.12) 
Jo 

where f(x) is the projection of the distribution Po(|^r|) onto the x axis 

f{x) = J dr V {\8r\,8t)8(x- \8r\) (5.13) 



which satisfies 



/oo 
dxx 2 f(x) = 2D8t . (5.14) 
-oo 



Equations (5.12) and (5.14) show that the mean inward displacement (Sx) x is of typical 
magnitude V DSt, in a layer of depth VD5t next to the boundary, and negligible elsewhere. 
The weight w of these inward displacements is clearly ~ DSt. We define 



w = lim 

L— >oo 



dx{5x) 



(5.15) 



We evaluate w by sustituting (5.12), then making a change of variables X = x' + x, 
X' = x 1 — x. The integral is written as the sum of two integrals, one over the domain 
X < L, \X'\ < X, which vanishes because of a symmetry, and another integral which 
involves only f(X') in the limit L — > oo: we find w = DSt, which implies that 



J dvF l {v){5r l ) = -DSt J ds l F l 



(5.16) 



for any vector field F, where dsi are the components of a vector element of the surface. 
When evaluating the integrals over (Sridrj) we can approximate these terms by 2D5t5ij, 
because the second term in (5.14) is significant only in a narrow layer of width \/D5t. The 
terms containing averages of dr 3 make no contribution at order St 2 . Retaining only the 
leading order terms, we find the following contribution for t ^ 0: 

C(t) = ^JdrJ dv'd r Mv)d rl A 3 {v')V{v,v'-t) 

-2^r J d Si j dv'A^d^A^V^v'-t) 

+ ^ JdsiJds'jAiWAjWVfas'it). (5.17) 
After integrating by parts, and adding the delta function contribution from t = 0, we find 

C(t) = D6(t) Jdr AiAi-^- J dr J dr' A { Aj d^. r ,V(r, r'; t) (5.18) 

Before Fourier transforming this expression to determine the absorption coefficient, we will 
introduce a convenient expression for the propagator: 



V(r, r'; t) = J] Xa(r) Xa W) exp(-Dk 2 J) 



(5.19) 



where the Xa( r ) are solutions of the Helmholz equation [A + fc„]xa = 0, satisfying the Neu- 
mann boundary condition riid ri Xa = 0. Using this result, we find the following expression 
for the absorption coefficient 



a(ui) = Klu z 



where 



K 



i) 2 + " 2 
ge 2 D E, 



B 2 



r2 • 

n 



(5.20) 



(5.21) 



Equation (5.20) is the main result of this section. As pointed out above, it consists of two 
terms. The first term is just the classical result derived at the beginning of this section. 
The second term is written as a sum over eigenmodes of the diffusion propagator. It 
incorporates boundary effects, as will be seen in the next section. 

5.3 Low-frequency limit 

We conclude this section with two remarks concerning the low frequency limit. First, 
we show that that the result (5.20) fulfils a condition discussed at the end of section 3, 
where we noted that the low-frequency limit of the absorption coefficient must be invariant 
under adding any gradient to the vector potential. To see this, consider the effect of the 
following transformation: 

Ai^Ai + d r .(p. (5.22) 

In the limit u — > the absorption coefficient is determined by the integral of the correlation 
function. The addition of the term d ri ip produces two new terms in the delta function 
contribution in (5.18), one an integral containing Aid ri ip, the other containing d ri fd ri (p. 
Corresponding additional terms appear in the double space integral in (5.18). Consider 
the first of these terms, linear in <9 r . ip. In order for these terms to make no contribution to 
the absorption coefficient, it is sufficient to show that the following two integrals are equal 

I = D J dt Jdrjdr' A, d r .<pdl. r ,V (5.23) 

/' = 2 J drAid ri ip. (5.24) 

We will use the result 

8(r-r') = J2xa(r)Xa(r') (5.25) 

a 

which follows from setting t = in (5.19). Using (5.20), we find 

/ = 2 5Zp" J dr J dr'Aidr^d^xM d r .x a (r') 

= 2 J dr J dr'Ai p(r') d r ,8{r - r') = I' . (5.26) 

This shows that terms involving Aid ri cp cancel and do not contribute to the absorption in 
the limit oj — > 0. Since this result applies for any vector field Ai, we can replace Ai by 
d ri ip and deduce immediately that the terms quadratic in d ri ip also cancel. 

Second, we comment on the form of the absorption coefficient in the limits oj ^> oj c and 
oj <C oj c . The vector potential can always be written as a sum of the curl of a divergenceless 
field, and a gradient 

\u) A = a + Vip = V A -0 + Vy? , V.^ = (5.27) 

with the field a chosen so that it is tangential to the boundary (n.a = 0). In the case 
u ^> o; c , the absorption coefficient is determined by the delta function contribution, and 
we have 

lim a(oj) = Ku 2 / dr (a 2 + V^ 2 ) (5.28) 



with K given by (5.21). In the case u <C u Cl on the other hand, the calculation we described 
above shows that the potential <p makes no contribution to the absorption coefficient, and 
that 



We showed in section 3 that if the vector potential is written in the form (5.27), then 
the potential ip is zero when the conductivity is isotropic and local. Inspection of (5.20) 
shows that the boundary contribution to the correlation function vanishes when ip vanishes. 
In more general cases, comparison of (5.28) and (5.29) shows that if the electric field is 
independent of frequency, the absorption coefficient is reduced, relative to its classical 
value, at frequencies below the Thouless frequnency u c . 




(5.29) 



6. Self-consistent electric field: diffusive case 



6.1 Non-local conductance 

We can deduce the non-local conductivity from the results of the previous section in 
two ways. We could use (4.6) as the definition of the conductivity, and evaluate it by 
setting A(r) = 5(r — R)ej in (5.5), so that 



e 2 



<7y(R,R';f) = ^ (8(K - v)6r i 8(^' - r>j) • (6.1) 

Alternatively, if we write the absorption coefficient in terms of the non-local conductivity 
in the form 

a(uj) = J dr J dr ' J drexp[ic^r]<Ty(r, r',u) Ei(r) Ej(r') (6.2) 

the kernel (Ty(r, r';t) is deduced from (5.18). By either route we find 

a(r, r'; t) = E k^(r - r')S(t) - Dd r 2 v P(r, r'; t)] (6.3) 

where E is the bulk conductivity, and V(r, r', t) is the propagator, given by (5.19). 

This form for the non-local conductivity was originally given in ref. [13] (and in the 
DC limit in ref. [12]). The argument in these earlier papers involves the diagrammatic 
analysis of disorder averaged perturbation theory, and in the case of [13] it appeals to a 
supersymmetric formalism. We believe that our derivation is more direct, and also more 
compelling. Our derivation considers the effect of the surface explicitly, whereas it is not 
clear from the diagrammatic analysis that there are not additional contributions which 
arise from integrating fields over the surface of the sample. Our derivation also deals 
explicitly with the fact that the trajectories are discontinuous, and we show explicitly how 
the correct evaluation is related to the Stratonovitch definition of the integral over the 
trajectory. 

6.2 Self-consistent solution 

The self-consistent electric field has to be chosen to satisfy (3.5). First we remark 
that in the case of diffusive electron motion, the term containing eo is negligible, and can 
be dropped. Estimating the magnitude of a by the bulk conductance E = ne 2 D, and 
noting that the bulk plasma frequency scale is cu 2 ~ Ne 2 /m e e 0: we see that this term is 
negligible provided to 3> cu s , where uj s is the elastic scattering rate. This is consistent with 
the assumptions that the electron motion is diffusive. 

We therefore wish to determine an electric field for which V.(aE) = 0; more explicitly 
we require 

E d n J dt exp(iwt) J dr' [Ei(r')8{r - r')5(t) - Dd 2 r yV{r, r'; t) Ej(r')] = . (6.4) 
Substituting for the propagator using (5.19), we find 



= E 



= E d ri Ei + E 

Of 

= Eq d Ti Ei — Eq 
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D 2 fc 



2lA 



D 2 k 



D 2 ki + u 2 



Xa (r) Jdr'd rj xUr')Ej(r') 
Xa {v)[jdY' Xa {v')d rl E 3 {v') 
-J dsixM Ei(s) 



(6.5) 



This equation has a solution where V. E = everywhere, with E tangential to the bound- 
ary. The classical solution for an isotropic local conductance satisfies these conditions. We 
conclude that, at least in the case where there is no static biasing magnetic field applied, 
the electric field configuration is independent of frequency. 

6.3 Calculation of the absorption coefficient 

We have shown that the electric field distribution is independent of frequency: we can 
therefore use the solution of the form (3.7), (3.8), with = 0, and with ip tangential to the 
boundary. Referring to (5.21) we observe that, after integrating by parts, the integrals are 
seen to vanish because the field A is divergenceless and is tangential at the boundary. It 
follows that the summation in (5.21) vanishes, and that the absorption coefficient is given 
by the classical expression (5.1), at all relevant frequencies. 

We conclude by describing a useful approach to calculating the field ip, satisfying 
Poisson's equation (3.9). The solution can be obtained from a Green's function G(r, r') 
satisfying V 2 G = —S(r — r'). A suitable Green's function is 



G(r,r')=£ 



Ur)Ur') 



U2 



(6.6) 



where the £ n ( r ) and /c 2 are eigenfunctions and eigenvalues of the Helmholtz equation, 
[V 2 + fc 2 ]£ n = 0, solved with the Dirichlet boundary condition £ n (r) = 0. The field ift(r) 
is then obtained by applying this Green's function to the source term iuB appearing in 
(3.9). In two dimensions, the absorption coefficient can then be written in terms of the £ n 
as follows 



a( U ) = J|/*|W| ! 



2El 



dr -0V 2 V 7 * = 



2Etf 



dr ip 



(6.7) 



Using the Green's function (6.6) to obtain ip, we have 



dr Cn(r) 



(6.8) 



This is a very general expression for the classical magnetic dipole absorption coefficient 
in particles with diffusive electron motion, expressed in terms of solutions of the two 
dimensional Helmholtz equation. As it stands, equation (6.8) is valid for two-dimensional 
particles. 

We note that the eigenfunctions £n( r ) obey Dirichlet boundary conditions and that 
the absorption coefficient for a given geometry can be significantly reduced by applying 
cuts orthogonal to the boundary: the main contribution to the sum in (6.8) comes from the 
ground state and the low-lying states and the corresponding eigenvalues are increased by 



applying cuts. This behaviour is expected since such cuts inhibit the flow of eddy currents 
which causes the absorption. 



6.4 Some examples for simple geometries 

In this section we summarize our results for specific geometries, namely discs, squares, 
and spheres (the result for the latter is well known [3], and is included to establish con- 
nections with earlier work). For two dimensional discs of radius a, we have ip = icoBor 2 /4. 
Using (3.10) and (5.1) we obtain 



a ^ = 16~^ — " ( } 



This result is easily shown to be consistent with (6.8), using the fact that the k n are defined 
by J (k n a) = and that J2 n ^n 4 = « 4 /32. For squares of sidelength a, we obtain 

. 32 Spc^a 4 ^ 1 1 
aUu) = ^ — > 9 — ^ ~ . (6.10) 

7T 6 c 2 m 2 n 2 m 2 + n 2 J 

m,n>0 

The sum can be evaluated numerically and gives 0.528 .... Finally, for spheres of radius 
a (6.8) is to be modified as follows 



2 .<". i .■• 2 



z u 2 r ^ i 

J -a _ r"n 



dr f n (r) 



2c 2 

n ■- - 

where the eigenvalues are defined by Jo(k n r±) = and r\ = a 2 — z 2 . We obtain 



(6.11) 



7T E w 2 a 5 

a(a;) = 15 c 2 ' ( 6 ^ 



Equation (6.11) reproduces the well-known absorption coefficient for metallic spheres [3]. 
In [3], this result is compared to the absorption coefficient for electric dipole absorption, 
which has a different size dependence, ~ a 3 . In two dimensions, on the other hand, the size 
dependence is the same for magnetic and electric dipole absorption. The latter coefficient 
was calculated in [9], 

, , 34 ela A u 2 
a el (u) = — -2—- ■ (6.12) 



9tt £ 



o 



7. Ballistic electron motion 



7.1 General remarks 

In the case where the electron motion is ballistic, the electric field must be determined 
by the non-local conductivity, and equations (3.5) and (3.6) must be solved to determine 
the electric field. Fortunately, in both the high and low frequency limits, there are consid- 
erations which simplify the discussion. 

In the low frequency limit, uj <C uj c , we showed in section 3 that only V A E is 
relevant, and that we can use any electric field for which the circulation is spatially uniform. 
It follows from (2.3)-(2.5) that the absorption coefficient is proportional to uj 2 at low 
frequencies. 

In the high frequency limit, the conductivity tensor (4.10) acts over a range R ~ v-p/u, 
and the conductivity becomes effectively local, with value S(cu) = \Ne 2 /m e uj. The non- 
local self consistency condition for the electric field then reduces to the same requirements 
as for the diffusive case: the electric field is tangential to the boundary, and is derived 
from a field ip which satisfies (3.9). 

^From (2.5), it is clear that the high frequency behaviour is determined by disconti- 
nuities in derivatives of the correlation function. The correlation function of the smooth 
perturbation f(t) = r.E for motion in a billiard has discontinuous derivatives due to the 
change of direction when the particle collides with the boundary. In the neighbourhood of 
a collision with the boundary at t = 0, the perturbation takes the form 

/(f) = (r' + r").E[(r' + r")t]6(-t) + (r' - r").E[(r' - r")t]Q(t) (7.1) 

where r' and r" are respectively tangential and normal components of the velocity at the 
instant before the collision, and 6(£) is the step function. Taylor expanding E(r) we find 
(with repeated indices summed over): 

/(f) ~f>E t + f'l Et {0(-f) - 6(f)} + (f# + r'/r'J) d T] E % t 

+{r' i r'; + r'lr' j )d rj E i t{Q{-t)-Q{t)} (7.2) 

For a general electric field, /(f) has discontinuities in t of magnitude 2f" E t on collision 
with the boundary, but for electric fields tangential to the boundary, the discontinuities 
are in the first derivative, and are of magnitude 2 {f^r" + r^'r'A d rj Ei. 

If fit) has discontinuities in its n th derivative, the Fourier transform of its correlation 
function decays as u; _2 ( n+1 ) as uj — > oo. The absorption coefficient is obtained from this 
Fourier transform by multiplying by a factor which contains uj 2 . In the case of a general 
field we therefore expect the absorption coefficient to approach a constant for uj ^> u c , 
whereas for a tangential field we expect that a(uj) ~ uj~ 2 for uj ^> uj c . We know that the 
field E(r, uj) approaches the tangential form as u — > oo, but we have no information about 
how rapidly this limit is appraoched. We can only say that the absorption coefficient must 
decrease for uj 3> uj c , and that it is unlikely to decrease faser than uj~ 2 . 

7.2 An example: the square billiard 

It is instructive to discuss an example: we consider the absorption coefficient for a 
square billiard with ballistic electron motion, with two different, frequency independent, 
choices for the electric field E(r), both satisfying (3.7), (3.9). First we calculate the 
absorption coefficient assuming that the field is tangential to the boundary, and then 
consider the case where it is circularly symmetric (which corresponds to taking the angular 



momentum operator L z as the perturbation. The results will illustrate the application of 
(2.5), and will verify two of the conclusions from the arguments above: we find that the 
absorption coefficients agree in the low frequency limit, and that at high frequencies the 
absorption scales as u~ 2 for the tangential field, but as uj° for the radially symmetric field. 
In terms of the perturbation 

Atf(r,p) = — p.A(r) (7.3) 

m e 

the absorption coefficient is given by 

a(u) = —g(E F )Re J dt e*"* (AH(r t , p t ) AH(r, p)) . (7.4) 

The perturbation (7.3) is determined by the choice of the vector potential A(r). We will 
first assume that A(r) is tangential to the boundary of the particle. Accordingly we take 
4> = in equation (3.7). The field ip(x,y) in equation (3.8) is determined from equation 
(3.9) which is most conveniently solved using the Green's function (6.6). For a square of 
side a we obtain for the perturbation 



Aff(r,p) = — B a 2 (-Y 

m e \7r/ 



E" <[ nn sin (mux /a) cos (niry/a) 



mn 



a mn(m 2 + n 2 ) 



\ mn cos (mux /a) sin (nny/a) 



mn 



a mn(m 2 + n 2 ) 



(7.5) 



Since motion in the square is integrable, the autocorrelation function of the perturbation 
in (7.4) is calculated as an average over tori 

AH{v u pt) AH(r, p)} = J -jLL d/x (i) AH(1, 6) AH(I, 6 + w(I)t) (7.6) 

where dfi(T) = g(Ep)~ 1 d[Ep — H(T)] averages over the tori. I and 6 are the action and 
angle variables characterizing the motion in the square, u?(I) are the respective frequencies. 
Equation (7.6) is easily evaluated [19]. The result is of the form 

8 m e e 2 uj 2 a 5 v F 

«M = — ^72 J W K (7-7) 

where f(z) is an energy-independent scaling function and u c = vp/a. For large frequencies, 
f(z) ~ z~ 4 and hence a(u) ~ u~ 2 . For small frequencies, on the other hand, one obtains 

„ o o c odd ^ 

_ 8 m e e z u A a°v F ^ _1 1 ,„ , 

a[UJ) ~^ c 2 h 2 ^ m 2 n 2 (m 2 +n 2 ) 3 / 2 ' 1 J 

m,n>0 

As remarked in section 7.1, the absorption coefficient is proportional to u 2 for small fre- 
quencies. 

In order to verify explicitely that the low-frequency absorption does not depend on the 
boundary conditions for the electric field, as discussed in sections 3.3 and 7.1, the above 
calculation can be repeated using a vector potential in the symmetric gauge 



A = -y(-y,x,0) 



(7.9) 



We note that this choice of the vector potential does note satisfy tangential boundary 
conditions. The corresponding perturbation is 



AH(v,p) = — p.A = ^L z . (7.10) 

m e 2m e 

For small frequencies we find again the result (7.8), thus verifying explicitly that the 
boundary conditions do not influence the low-frequency absorption. For high frequencies, 
on the other hand, we find a(u) ~ cu°, as predicted in the previous section. 
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Figure captions 



Figure 1: Illustrating the images used in discussing expression (4.11) for a square billiard. 
The sum over all paths from r to r' can be represented as a sum over straight lines k of 
length 1^ connecting r with the image points of r'. 

Figure 2: Illustrating the vectors and coordinate system used in the discussion of the 
construction of the short time propagator. 



